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Abstract 

We review the phenomenology of coupled magnetic and electric order parameters for systems in 
which ferroelectric and incommensurate magnetic order occur simultaneously. We discuss the role 
that such materials might play in fabricating novel magnetoelectric devices. Then we briefly review 
the mean-field description of ferroelectricity and modulated magnetic ordering as a preliminary to 
analyzing the symmetry of the interaction between the spontaneous polarization and the order 
parameters describing long-range modulated magnetic ordering. As illustration we show how this 
formulation provides a phenomenological explanation for the observed phase transitions in Ni3V208 
and TbMnOs in which ferroelectric and magnetic order parameters simultaneously become nonzero 
at a single phase transition. In addition, this approach explains the fact that the spontaneous 
polarization only appears along a specific crystallographic direction. We analyze the symmetry of 
the strain dependence of the exchange tensor and show that it is consistent with the macroscopic 
symmetry analysis. We conclude with a brief discussion of how our approach might be relevant in 
understanding other systems with coupled magnetic and ferroelectric order, and more importantly, 
how these principles relate to the search for materials with larger magnetoelectric couplings at 
room temperature. 

PACS numbers: 



I. INTRODUCTION 



The interactions between long-range magnetic order and long-range ferroelectric order 
have been studied in depth since the first experimental confirmation of the magnetoelectric 
effect in the late 1950s.i2i^ We note the existence of several reviews^i^ and monographs^ 
which give a general overview of the subject. 

Of particular interest for this review are those materials which exhibit a combined 
magnetic and ferroelectric transition. Perhaps the best known of these is Ni-I boracite 
(NisByOial) which shows coupled ferromagnetic, ferroelectric, and ferroelastric properties 
at a single phase transition at T = 61.5 K.i^ The multiferroic behavior in this boracite arises 
from the fact that the magnetic transition is connected to a structural distortion, which in 
turn allows the development of ferroelectric order.^^ This transition can be understood in 
terms of a phenomenological Landau theory which couples the ferromagnetic, ferroelectric, 
and ferroelastic order parameters to a primary antiferromagnetic order parameter.— The 
strong coupling between magnetic and ferroelectric order parameters in systems having a 
simultaneous phase transition is demonstrated by the observation that in NisByOial it is 
possible to reverse the direction of the spontaneous polarization by applying an external 
magnetic field perpendicular to the direction of magnetization.— 

Cr2Be04 also develops magnetic and ferroelectric order at a single phase transition.— 
Below T=28 K, Cr2Be04 orders antiferromagnetically into a state with spiral spin struc- 
ture, and this antiferromagnetic state shows an extremely small spontaneous polarization 
(approximately one million times smaller than that of BaTiOa). The coupling between mag- 
netic and ferroelectric order is expressed by a model proposing a mechanism in which the 
electric polarization is induced solely by the antiferromagnetic order.^ A similar model for 
magnetically-induced ferroelectric order will be discussed in detail in the following sections. 

While the magnetic and ferroelectric transition temperatures for BaMnF4 are widely sep- 
arated, this system is useful in illustrating the importance of symmetry considerations in 
determining magnetoelectric properties. Pyroelectric BaMnF4 orders antiferromagnetically 
when cooled below T = 26 K, and there is a dielectric anomaly at this magnetic transi- 
tion temperature.— This decrease in dielectric constant below T/v varies like the square of 
the sublattice magnetization, and clearly indicates a coupling between the magnetic and 
ferroelectric properties of the sample. This interaction between magnetic and ferroelectric 
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order is attributed to a magnetoelectric coupling which causes a polarization induced spin 
canting.— Substituting l%Co for Mn changes the magnetic symmetry group of the com- 
pound to one which precludes this magnetoelectric coupling,^- and in turn eliminates the 
dielectric anomaly at the magnetic ordering temperature. As this system illustrates, in 
order to understand magnetoelectric couplings in multiferroic systems it is crucial to have 
complete information about the magnetic and structural symmetries of the system. 

Until quite recently, the theoretical and experimental studies have focussed on ferroelec- 
tricity in systems with simple ferromagnetic or antiferromagnetic orderi^iii (with studies on 
Cr2Be04 being the notable exception). These systems are tractable from a theoretical stand- 
point, and allow a comparison to be made between experimental results and straightforward 
models based on magnetic space groups. However, limiting the scope of investigation to 
systems with ferromagnetic or antiferromagnetic order neglects a large class of materials 
which have more complex magnetic structures. Here we will not consider systems (several 
of which are listed in Table I of Ref. 0) which are ferroelectric at high temperature and then 
have a lower temperature phase transition at which magnetic ordering takes place.-^ Instead, 
in this brief review article we will focus on the more recent studies in which ferroelectricity 
appears simultaneously (in a single combined phase transition) with long-range sinusoidally 
modulated magnetic order,"lSJi which we will refer to generically as "incommensurate"— 
magnetic order. Accordingly, we will briefly summarize the experimental situation for the 



systems TbMnOa (TMO)™ and NigVaO^lNVO) 



the symmetry analysis developed in Refs. I2C 
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^ Then we will describe in detail 



23l to understand the phenomenol- 



ogy of these systems. We believe that this theoretical approach is simple enough that it 
can easily be applied to the ever increasing number of systems like NVO or TMO in which 
ferroelectricity is induced by incommensurate magnetic long-range order. 
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To illustrate this phenomenon, we show, in Fig. ^ some intriguing data from Ref. 
showing that the spontaneous polarization P depends strongly on the applied magnetic field 
H. At first glance this data seems to have no obvious explanation. However, when viewed 
in combination with the magnetic phase diagram (see Fig. IHl below) we will see that this 
data indicates that the spontaneous polarization is nonzero only in the magnetic phase we 
will call the "low temperature incommensurate phase." The hysteresis is a consequence of 
passing through a first order phase boundary between this phase and an antiferromagnetic 
phase in which a spontaneous polarization is not allowed. Thus the dramatic dependence of 
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FIG. 1: Adapted from Ref. Spontaneous (i. e. in zero applied electric field) polarization 
P versus applied uniform magnetic field H applied along different crystallographic directions at 
selected temperatures for NVO. The arrows indicate the directions of increasing and decreasing 
magnetic field. 



polarization on magnetic field has a simple explanation: ferroelectric order appears only in 
one specific magnetic phase whose existence depends in the value of the magnetic field. This 
strong coupling between magnetic and ferroelectric order is potentially important for device 
applications, as we will discuss in the following section. From a basic physics standpoint, 
these systems which exhibit a coupling between the ferroelectric moment (a polar vector) 
and the magnetic moment (an axial vector) are very interesting. (As we will see, such 
systems have order parameters whose response to both electric and magnetic fields becomes 
large especially near a phase transition.) A complete understanding of this coupling from 
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a microscopic theory is not yet available. Here we will show that the Landau expansion 
explains the observed phenomenology of this interaction and that these results follow from 
the microscopic symmetry of the strain dependence of the exchange tensor. This explanation 
will serve as a guide to constructing a fully microscopic theory of magnetoelectric coupling. 

Briefly, this review is organized as follows. In Sec. II we discuss some general types of 
applications in which the magnetoelectric coupling may be exploited to develop new types of 
devices. It should be emphasized that these applications are speculative, and are intended 
to illustrate the types of new devices that could be developed using these new materials. In 
Sec. Ill we review the Landau description of ferroelectricity. In Sec. IV we give a simplified 
theoretical analysis of incommensurate magnetic ordering and in Sec. V we discuss how 
Landau theory leads to a symmetry-based description of incommensurate magnetic order- 
ing. It is our aim to demystify the use of representation theory for the determination of 
magnetic structure by diffraction techniques. Understanding these incommensurate mag- 
netic structures is crucial to developing a model for the coupling between magnetic and 
ferroelectric order in these systems. In Sec. VI we use the results of Sec. V to analyze how 
symmetry restricts the form of the coupling between electric and magnetic order parameters 
and thereby explain the simultaneous appearance of these two kinds of order parameters in a 
single phase transition. The construction of this interaction is greatly simplified by the fact 
that it involves an expansion in powers of the order parameters relative to the paramagnetic 
paraelectric phase. Thus the interactions have to satisfy the invariances of the disordered 
paramagnetic/paraelectric phase2^»2^ and we do not need to broach the more complicated 
question of analyzing the symmetry of interaction within an ordered phase. In Sec. VII 
we analyze the symmetry of the strain dependence of the exchange tensor and show that 
it leads to results identical to those of Landau theory. Finally in Sec. VIII we summarize 
the main points of this review and speculate on some future directions of research. We will 
discuss how our results on ferroelectric order in incommensurate magnets may offer guidance 
in searching for new magnetoelectric materials. 



II. DEVICE APPLICATIONS 



The development of devices incorporating both charge and spin degrees of freedom, often 
referred to as spintronics, has already led to significant technological breakthroughs.-^ Mag- 
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netic sensors based on giant magnetoresistance (GMR) are widely used as the read heads in 
modern hard drives, and magnetic random access memory also relies strongly on couplings 
between charge and spin. Additionally, there are a wide range of proposals for devices based 
on controlling the spin degree of freedom in ferromagnetic semiconductors, including spin 
valves and qubits for quantum computing. Much of the research on materials in which charge 
and spin are coupled have focussed on metallic and semiconducting systems. However, di- 
electric materials exhibiting couplings between electric polarization and magnetization may 
also play an important role in developing the next generation of spintronic devices. 

Magnetoelectrics are systems in which either applying an external magnetic field produces 
an electric polarization or applying an external electric field produces a magnetization. This 
type of coupling between charge and spin was postulated by Pierre Curie at the end of the 
19th century, ^"^ but not observed experimentally until the late 1950s.^'^ Materials in which 
two or more of ferroelectric, ferromagnetic, and ferroelastic order coexist are referred to 
as multiferroics. This strict definition of multiferroics is often relaxed to include systems 
which exhibit combinations of any type of long range magnetic, ferroelectric, or ferroelastic 
order. This review will concentrate specifically on magnetoelectric multiferroics, where the 
coexistence of long range magnetic and long range dielectric order leads to a pronounced 
couplings between the charge and spin degrees of freedom in these systems. 

We consider two classes of devices based on magnetoelectric multiferroics. The first class 
of devices depend on the magnetoelectric effect — the induction of a magnetization (polariza- 
tion) by an applied electric (magnetic) field. Using the magnetoelectric effect, it is possible 
to design a range of devices from sensors to transducers to actuators, coupling magnetic 
and electric properties. The second class of devices exploits the fact that these materials 
have simultaneously appearing long range magnetic and ferroelectric order. The underly- 
ing assumption is that multiferroics exhibit both charge and spin ordering, and due to the 
coupling between the two, both magnetic and ferroelectric order will be strongly affected by 
either magnetic or electric fields. Strictly speaking, only magnetic field control of the elec- 
tric polarization has been demonstrated for the multiferroic materials with incommensurate 
magnetic structures discussed in this review, but magnetic phase control by an electric field 
has been demonstrated in other multiferroic materials. This coupling between long range 
electric and magnetic order leads to new functionalities which can be exploited for designing 
new types of spintronic devices. 
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The investigation of magnetoelectric devices is an active area of research. Prototype de- 
vices fabricated using piezoelectric-magnetostrictive composite materials to produce magne- 
toelectric coupling have already been tested,—"^ and there are a range of proposals for other 
magnetoelectric devices. These include utilizing magnetoelectric materials as the pinning 
layer in GMR devices^- ^ for low frequency wireless power applications,™ and for developing 
tunable dielectric materials.!^ One key feature of magnetoelectric materials is that they al- 
low the design of devices controlled magnetically or electrically, as desired. Controlling the 
magnetic properties of materials using an electric field offers significant benefits in designing 
new devices. Using current-based methods to switch magnetic devices is relatively slow, and 
power-intensive. Voltage control of the magnetic properties is expected to offer significantly 
faster switching (thin film ferroelectrics can show switching times of less than 200 pa^) 
in a low-power device. Magnetoelectric materials offer the potential for fabricating highly 
tunable, fast switching, low- loss/low-power devices having very small form factors, which 
would be suitable for a wide range of commercial and industrial applications. 

The materials property most relevant in determining the suitability of a compound for ap- 
plications in magnetoelectric devices is the magnitude of the magnetoelectric susceptibility, 
Xme- For homogeneous materials, Xme satisfies the bound. 



where xe and xm are the electric and magnetic susceptibilities of the system respectively.- 
Therefore, in order to maximize the magnitude of the magnetoelectric coupling, one should 
attempt to maximize the magnitudes of both xm and xe- Since ferroelectrics typically have 
large values of xe and ferromagnets typically have large values of xm-, multiferroics are 
expected to have large values of Xme- Furthermore, since susceptibilities are largest at the 
ordering transition, systems developing magnetic and ferroelectric order at the same temper- 
ature should show exceptionally large magnetoelectric couplings. This has been confirmed 
for the intrinsic multiferroic Bao.5Sri.5Zn2Fei2022, which has the largest magnetoelectric 
coefficient of any single-phase material identified to date.»24 Understanding the microscopic 
origins of the magnetoelectric coupling in these multiferroic systems will have important 
ramifications for developing novel magnetoelectronic devices. 

Beyond simply exhibiting very large magnetoelectric couplings, intrinsic multiferroics 
also have both long range magnetic order and long range ferroelectric order. The coupling 




(1) 



7 



between magnetization and polarization offer new possibilities for designing devices. The 
ability to control the magnetic or ferroelectric state of a system using either a magnetic field 
or an electric field would offer the ability to develop multifunctional memory elements, for 
example, ferroelectric memory which can be written to using magnetic fields. We will discuss 
two proposals for new technologies which explicitly utilize the ferroelectric and magnetic 
characters of magnetoelectric multiferroics. It should be emphasized that this discussion is 
meant only to illustrate some of the potential applications arising from the incorporation 
of multiferroic materials into new devices. More investigation on the specific properties of 
these multiferroics is required before proof-of-principle devices could be designed based on 
these speculations. 

As the bit density of modern hard drives increases, the characteristic size of the mag- 
netic structures used to store the information is decreasing. As the physical size of the bit 
is reduced, the anisotropy energy decreases, and the magnetic moment can begin to ther- 
mally fiuctuate. Controlling these thermal fiuctuations is necessary to ensure the long-term 
stability of stored information in ultra-dense magnetic recording material. For long term 
magnetic storage (5-|- years), the ratio of the energy barrier against these thermal fluctua- 
tions to ksT should be large, roughly 50. In current devices, this is often accomplished by 
using materials with very large magnetic anisotropy energies or by exploiting the anisotropic 
difference between FM and AFM layers. One possible apphcation for multiferroic materials 
is to use the coupling between ferroelectric and magnetic order in these systems to stabilize 
the magnetic moment against thermal fluctuations in nanoscale magnets. 

In many magnetoelectric multiferroics there is a strong coupling between the ferroelectric 
and magnetic order parameters. In such systems, fixing the polarization (magnetization) 
direction will fix the axis of the magnetization (polarization). This coupling is observed in 
measurements showing that the sign of the magnetically induced polarization is indepen- 
dent of the sign of the applied magnetic field, although the development of ferroelectric 
order depends strongly on the magnetic field axis. In such multiferroics, fixing the electric 
polarization would also fix the magnetization axis. This ferroelectrically induced magnetic 
anisotropy would inhibit thermally activated switching of the magnetic moments by signifi- 
cantly increasing the magnitude of the energy barrier to magnetization reversal. This could 
be accomplished, for example, by assembling multiferroic nanoparticles on a ferroelectric 
substrate. In this geometry, the very large ferroelectric anisotropy energy would provide a 
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tunable barrier against thermal fluctuations of the magnetic moment as well. 

Multiferroics may also have important applications in developing magnetic fleld sensors. 
There are a range of proposals for incorporating magnetoelectric materials in exceptionally 
sensitive magnetic fleld detectors. Even relatively small external magnetic flelds will produce 
a voltage change in materials with very large magnetoelectric couplings. Since it is often 
better to measure small voltages at zero applied current rather than small magnetizations 
or small changes in resistivity, magnetoelectric materials offer the potential for developing 
greatly improved magnetic fleld sensors. Because multiferroics exhibiting simultaneous mag- 
netic and ferroelectric transitions offer exceptionally large magnetoelectric couplings, these 
materials are particularly interesting in the context of improved sensors. Figure |21 shows 
a schematic for such a device.— The magnetization produces a spontaneous polarization 
directed perpendicular to the plane of the sensor. This magnetically induced voltage can 
be measured to a high degree of accuracy, either directly, or by measuring the dielectric 
response of the compound. This device could also be conflgured to extract energy from an 
alternating magnetic fleld — the magnetically induced alternating voltage could be used as a 
supply for very low power applications.— 

Beyond simply being used as a passive magnetic fleld sensor, the device illustrated in 
Fig. 121 could also be conflgured as a voltage biased magnetic memory element. One of the 
difficulties facing current magnetic random access memory (MRAM) devices lies in producing 
sufficiently strong magnetic flelds to cause a moment reversal in the memory element, but 
also sufficiently localized to affect only one specific element. While identifying multiferroic 
materials in which applying a voltage could reverse the direction of the magnetization would 
certainly be beneflcial for developing MRAM devices, a more modest type of voltage-assisted 
magnetization reversal could also be signiflcant. As will be discussed in the following section, 
ferroelectric order can be promoted or suppressed by the application of an external magnetic 
fleld in many multiferroic materials. We expect that in these materials, applying an electric 
fleld could then suppress or promote magnetic ordering. In such a system, the coercivity of 
the magnetic memory element could be tuned by applying an electric fleld. Consider the 
multiferroic memory element in a ferromagnetic state, which can be suppressed by applying 
a sufficiently large voltage. In the absence of an electric field, the coercivity of the memory 
element is large, so the magnetization is unaffected by stray magnetic fields. In order to 
reverse the magnetization direction, a bias voltage is applied to the multiferroic, bringing 
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FIG. 2: Schematic illustration of a device to measure magnetic fields by using the induced polar- 
ization. The middle layer (white) is a multiferroic material with strong magnetoelectric coupOlings 
and the outer layers (gray) are ferromagnetic metals. In this example the polarization is induced 
in a direction perpendicular to the applied magnetic field. 



the system closer to the magnetic transition, reducing the magnitude of the coercive field. 
In this state, the magnetization can be reversed by a relatively small external magnetic field, 
smaller than the coercive field of the unbiased multiferroic. When the voltage is removed, 
the new magnetization will be stable. This type of voltage-assisted magnetization reversal 
could be used to produce arrays of magnetic memory elements which could be switched by 
the same external magnetic field. Only those elements which have a bias voltage applied will 
have a sufficiently small coercivity to be switched by the magnetic field. This technique may 
offer advantages over transitional MRAM devices, such as a smaller sensitivity to stray fields 
(allowing higher bit density) and potentially faster switching times. This is schematically 
illustrated in Fig. |21 
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FIG. 3: Voltage-assisted MRAM. Left: Schematic diagram showing how the coercive magnetic field 
is expected to vary with voltage. Right: With the magnetization of the bits originally directed to 
the right, a small (50 Oe) writing field is applied. At zero voltage, this is smaller than the coercive 
field, so that the magnetization does not switch. However, applying 5V across the multiferroic layer 
reduces the coercive field to almost zero, so that this small magnetic field is sufficient to reverse 
the direction of magnetization. 



III. FERROELECTRICITY 



We start by making a few observations concerning the symmetry properties of ferro- 
electric systems for which magnetic ordering plays no role. In the most common scenario, 
ferroelectrics exhibit a high-temperature phase having spatial inversion symmetry which 
prevents the existence of a vector order parameter. Then, as the temperature is reduced 
through a critical value, Tp, a lattice instability develops in which inversion symmetry is 
broken cooperatively via a continuous phase transition at which a spontaneous polarization 
appears. Within a Landau theory this transition is described by a free energy of the form 

F = Ixe'P' + 0{P)' = ^a{T - Tf)P' + 0{P)' . (2) 

At the transition the fact that the quadratic term in P becomes unstable (negative) reflects 
the divergence in the electric susceptibility at the ferroelectric transition. This instability 
is sometimes traced to a soft phonon, but whatever the mechanism, the appearance of 
ferroelectricity represents a broken symmetry. Conversely, as will become relevant in the 
following, ferroelectricity can only occur if the symmetry is broken to permit the ordering 
of the polarization vector. We will use this criterion to determine which types of magnetic 
order can possibly induce ferroelectric order. If one takes the quartic terms in Eq. (0) to 
be of the form a[P^]^ (with a > for stability), then minimization of F with respect to P 
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shows that for T < Tp one has P ~ [2a(Tp — T)]^/^, which is expected to hold as long as 
Tp — T is not so large that sixth and higher order terms in F are important. Mean field 
theory ignores spatial correlations which lead to modifications of critical exponents, but the 
scope of this review does not permit consideration of such corrections.— 

As the temperature is further lowered it is possible for this ferroelectric system to develop 
long-range magnetic order.— In this case, one does not expect significant interaction between 
electrical and magnetic properties because the two phenomena are essentially independent 
of one another. In these systems, the spontaneous polarization will depend only weakly on 
the applied magnetic field. In this scenario, it is well known'^^ that one can expect anomalies 
in the dielectric response of the system when the ferroelectric develops (independently) long- 
range magnetic order. This review is not concerned with such an "accidental" superposition 
of electric and magnetic properties. Instead we focus our attention on the situation when 
the appearance of long-range magnetic order induces ferroelectricity. Furthermore, we will 
consider an interesting subclass in which the long-range magnetic order is modulated with 
an apparently incommensurate wavevector. We will develop a Landau theory for this com- 
bined phase transition in which the fact that the wavevector does not have high-symmetry 
(and is thus neither ferromagnetic or antiferromagnetic) is crucial to our analysis. Thus the 
development here can not be obtained by a trivial extension of theories applicable to ferro- or 
antiferromagnetic ferroelectrics. A simplifying feature of this formulation is that it is based 
on an expansion of the free energy in powers of the various order parameters relative to the 
paramagnetic phase. Accordingly, each term in this expansion has to have the full symme- 
try of the disordered phase.— In contrast, it is less straightforward to analyze whether 
or not the symmetry of a magnetically ordered phase permits an induced ferroelectric or- 
der. Also, the Landau formulation correctly predicts which components of the spontaneous 
polarization vector are induced by the magnetic ordering. In addition, the Landau expan- 
sion indicates that the spontaneous polarization is, crudely speaking, proportional to the 
emerging magnetic order parameter. 
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IV. TOY MODELS FOR INCOMMENSURATE MAGNETISM 



A. Review of Mean Field Theory 

In this section we review the description and phenomenology of incommensurate magnets, 
because the characterization of their symmetry is essential to understanding the coupling 
between magnetic and electric long range order. 

For the purposes of this review it suffices to consider the description of incommensurate 
magnets within mean field theory. For a system consisting of quantum spins of magnitude 
5* on each site, we write the trial free energy is 

F = U -TS = TT[pn + kTp\np] , (3) 

where 7i is the Hamiltonian, T the temperature, U the internal energy, S the entropy, and 
the actual free energy is the minimum of F with respect to the choice of p subject to the 
conditions that p is Hermitian with unit trace. Within mean field theory we take the density 
matrix to be the product of independent single particle density matrices p{i) for each site i: 

P^IIpW- (4) 

i 

This approximation corresponds to the intuitive idea that when correlations between spins 
are neglected, each spin reacts to the mean field of its neighbors. 

In Eq. (jni) the trace of pTi. gives the internal energy U and that of —kplnp gives the 
entropy S. In the absence of anisotropy it suffices to set 

P(0 = ^^[T + ccT{i) . S.] , (5) 

where X is the unit matrix of dimension (25" + 1), c is a constant of order unity, chosen 
to make Eq. (6) true, and Sj is the vector spin operator for site i [Here Sj is a (25* + 1) 
dimensional matrix] . The free energy is then minimized with respect to the trial parameters 
cr(i), which physically are identified as the average spin vectors: 

(S(z)) ^ Tr[p(z)S(z)] = a{{) . (6) 

Thus cr{i) is the vector order parameter at the ith lattice site. In this formulation the 
internal energy is quadratic in the order parameter cr, whereas the entropic term involves 
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both quadratic and higher powers of the order parameter. As we shall see, even without 
explicit calculations much information can be inferred from the symmetry of the trial free 
energy as a function of the order parameter (s). 

As mentioned in the introduction, we will focus our attention on systems which display 
incommensurately modulated magnetic long range order. We refer the reader to a compre- 
hensive survey of such systems by Nagamiya.— Here we give a simplified review. To char- 
acterize an incommensurate state we consider a toy model consisting of a one dimensional 
system with isotropic antiferromagnetic exchange interactions Ji and J2 between nearest 
and next-nearest neighbors, respectively. If J2 is antiferromagnetic and large enough, these 
two interactions compete and produce an incommensurate spin structure. Thus we are led 
to consider the Hamiltonian 

'H = '^^Sn ■ [JlSn+i + J2^n+2\ , (7) 

n 

with J2 > 0. The corresponding trial free energy is 

F = \dkT + J^a{i) ■ a{i + n) + 0{ct^) , (8) 

^ i n=l,2 

where the entropic term is scaled by a constant of order unity, d. 
B. Wavevector Selection 

It is instructive to write the free energy per spin, /, in terms of Fourier variables, cr(g) = 
(1/A^) I]j e*'''^''<T(z), where is the total number of spins as 

f ^ F/N=W X{q)-'^{q) ' ^{-q) + 0{a^) , (9) 

where x{q)^^ = dkT + Ji cos(ga) + J2Cos(2ga) is the wavevector-dependent susceptibility. 
At high temperature (when kT \Ji\ and kT ^ |^2|); x{q)~^ is positive for all q and the 
free energy is minimized by setting all the order parameters cr{q) to zero. In Fig. E]we show 
x{q)~^ a function of wavevector q for a sequence of temperatures. As the temperature is 
lowered through a critical value Tc, x{q)~^ becomes zero for the wavevector q = q^ which 
minimizes x{q)~^'- 

cos(goa) = -Ji/(4J2) • (10) 
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FIG. 4: x{q) ^ at a sequence of temperatures, T = Tc + 0.3n for n = 0, 1,2,3,4,5 with Ji = 1.0 
and J1/J2 = 2.56 (as is appropriate for NVO). Here Tc = 0.711. 

This determination of the value of go is called wavevector selection. As the temperature 
is reduced through the paramagnetic phase becomes unstable against the formation of 
long range order at the selected wavevector go- That is, for T < Tc the order parameter 
cr(go) assumes a nonzero value determined by the (negative) quadratic terms in combination 
with the (positive) terms of order cr^, so that |cr(go)| ~ (T^ — T)^^"^. Once order develops 
at one wavevector, the terms of order cr^ prevent order developing at other wavevectors. 
This scenario is realistic for a three dimensional system (for which long-range order is not 
destroyed by thermal fluctuations). The eigenvector associated with the eigenvalue of the 
quadratic form which passes through zero is called the critical eigenvector. The critical 
eigenvector contains the form factor of the ordering, i. e. it completely describes the 
pattern of spin ordering within a unit cell. In this simple model there is only one spin per 
unit cell, so the eigenvector specifies the direction i. e. the component which condenses. 
(This concept will be better illustrated when we consider real systems which often have 
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more than one magnetic site per unit cell.) In the present case when there is no anisotropy, 
the spin structure when cr(g) becomes nonzero for g = go is a modulated one in which the 
x-component of spin has a complex amplitude, A^, so that 
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and similarly for the other spin components. If these complex amplitudes all have the 
same phase,-^ then the spin is linearly polarized with an amplitude which varies sinusoidally 
with position. If the complex amplitudes do not have the same phase, then the spin structure 
will be a helix, a spiral, or a fan, etc. 

C. Effects of Anisotropy 

This toy model will not accurately capture the behavior of real magnetic systems because 
we have not yet included any anisotropy. In the presence of single-ion easy axis anisotropy, 
the trial free energy at quadratic order assumes the form 



where K is an anisotropy energy which favors alignment of spins along the easy axis, here 
the X-axis, and / denotes the free energy per spin. In this case, the instability (at which 
long-range order first appears) is one in which the spins are confined to the easy axis and 
have a sinusoidally varying amplitude. This type of ordered phase will be referred to as 
the high-temperature incommensurate (HTI) phase and the associated critical temperature 
will be denoted Thti- If the anisotropy is not too large, then, as the temperature is further 
reduced, the fourth order terms in the free energy (which we have so far neglected) become 
important. One effect of these terms can be visualized as incorporating the constraint of 
"fixed length." In the HTI phase the spins have a length which varies sinusoidally with 
position. However, in the ground state, we expect each spin to have its maximum length 
S but to be oriented in a direction to optimize the energy. Thus, in the extreme limit 
of zero temperature, the constraint of fixed spin length is fully enforced. Although the 
constraint is less fully realized at higher temperature, the qualitative effect is clear: when 
the temperature is sufficiently reduced, one has a continuous phase transition into a phase 
we refer to as the low-temperature incommensurate (LTI) phase. In this phase the spins 



(12) 




16 



Q n l~ I I I I I I I I I i I 1 I I I I I I I I I I I I 

0.00 0.25 0.50 0.75 1.00 1.25 
K/J^ 

FIG. 5: Zero temperature phase diagram for the J1-J2 {S = 1) model with easy-axis anisotropy 
scaled by The points represent numerical implementation of mean-field theory except for K = 
where analytic results are used. Here "HTI" denotes a longitudinally polarized incommensurate 
phase and "LTI" an elliptically polarized incommensurate phase. In both phases the modulation 
vector is given by Eq. (jlUj) . "AF" denotes a two-sublattice collinear antiferromagnetic phase. 
For large anisotropy this model reduces to the anisotropic nearest next-nearest neighbor Ising 
(ANNNI) model4a The dashed hue is dr awn for a value of the anisotropy energy which reproduces 
the evolution of magnetic phases in NVO as a function of T for H = 0. 

develop transverse order (in addition to the preexisting longitudinal order along the easy 
axis) to more nearly achieve fixed spin length. If the easy axis anisotropy is small, the 
range of temperature over which the HTI phase is stable is also small. The phase diagram 
of such a model as a function of anisotropy energy K and temperature is shown in Fig. 

We will mainly be concerned with the two incommensurate phases, the longitudinally 
modulated HTI phase and the elliptically polarized low-temperature incommensurate LTI 
phase. Although the details of the unit cell complicate the picture, the phenomenology of 
the HTI and LTI phases are usually roughly similar to that of the simplified case discussed 
here. In Fig. IHlwe show the experimentally determined phase diagrams of NVO and TMO 
as a function of applied magnetic field H and temperature T, with the HTI and LTI phases 
labelled. 
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FIG. 6: Phase diagram for NVO (left) and TMO (right) as a function of H and T. The HTI and 
LTI phases will be discussed in more detail below, but correspond roughly to the scenario of our 
toy model. The interpretation of the i?-dependence of the phase boundaries for NVO is given in 
Ref. l22l . For NVO the C and C phases are commensurate antiferromagnetic phases with a weak 



ferromagnetic moment. The NVO phase diagram is adapted from Ref. 
diagram is adapted from Ref. l4ll . 
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and the TMO phase 



D. Wavevector Locking 

From Eq. pO|) it would seem that the wavevector g is a continuous and smooth function 
of J2I Ji- Although our toy model does not give any simple explanation for the observed 
temperature dependence of g, a more complete analysis [as in Refs. and Q] leads to a 
small dependence on temperature which, like the dependence on J2/ Ji-, might be thought to 
be smooth and continuous. However, there are terms which favor commensurate values of 
q. These terms in the free energy must conserve wavevector, but only to within a reciprocal 
lattice vector G (which for our one- dimensional toy model can assume the values = 
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{277171 /a), where a is the nearest-neighbor separation). Thus one has the so-called Umklapp 
terms such as 6f = wa{qY, when 4g = Gm- More generally the Umklapp terms give a 
contribution to the free energy of the form 

= ^W^^n(T{qf''5nq-2m-K/a , (13) 

m,n 

where the coefficient Wm is of order unity and 5x is unity if x = and is zero otherwise. 
(Within the present formulation these terms come from expanding Trplnp to higher than 
quadratic order in the order parameter.) The effect of these Umklapp terms is to cause 
the wavevector to "lock" onto a commensurate value q = 2(m/n)7r/a as J2/J1 is varied. 
Since a{q) is smaller than one, especially near the ordering transition, these terms become 
much less important as the integer denominator 71 increases. Thus the effect of the Umklapp 
terms is that the variation of g as a function of a control parameter (such as the temperature) 
becomes a so-called Devil's staircase, which may either be complete or (if a is small enough) 
incomplete, as shown in Fig. [7| In the systems we will discuss here, there is no clear evidence 
of a Devil's staircase as a function of temperature. Accordingly, we find it convenient to 
imagine that q is incommensurate, and does not get stuck on commensurate values by 
Umklapp terms. Even if this is not strictly true, the difference in properties between an 
incommensurate system and a commensurate system with a large integer denominator n is 
experimentally irrelevant for the large values of n (n ~ 50), for the systems we will discuss. 
Accordingly, we will refer to the systems as "incommensurate" even though this may not be 
strictly accurate. 

In principle, the symmetry of real systems is usually such that anisotropy also occurs in 
the exchange interaction, in which case the trial free energy assumes the form 

/ = E x(g);V„(g) ■ ^„(-g) + 0(<T^) , (14) 

q a=x,y,z 

where x(g)~^ = dkT + Jiq, cos(ga) -|- J2Q, cos(2ga). If J„q, is isotropic (i.e. if it does not 
depend on a), then the wavevector selected for the ordering of the a component of spin also 
will not depend on a. However, in principle Jna depends weakly on a, and therefore the 
selected wavevector go will also depend weakly on a and the ordering will involve <7x{qx), 
ay{qy), and crz{qz)- Thus in the LTI phase it is possible that the two components of spin 
might have slightly different wavevectors, which we denote gnxi and ^lti- But as with the 
Umklapp contributions, there will be quartic terms in the free energy (in this formulation 
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FIG. 7: The Devil's staircase (DS). Left: an incomplete DS for the dependence of the magnetization 
of TbNi2Ge2 as a function of magnetic field4^ Right: Schematic example of a complete DS function 
Y = f{X). In a complete DS the function is nonconstant on a set of measure zero. For an 
incomplete DS the function is nonconstant on a set of nonzero measure. 

coming from the entropic terms) which favor locking the two wavevectors to be equal. These 
terms can be of the form 



where ciHTie*'^"'^^ = cr-iYniown) (cltiC*"'^™ = ctltiIq'lti)) is an order parameter characteriz- 
ing the appearance of the HTI (LTI) phase, and for simplicity we have assumed that the 
constant a is real-valued. This interaction only satisfies wavevector conservation if the two 
wavevectors are exactly equal. If, in the absence of this term, the two wavevectors are suffi- 
ciently close to one another, then this locking energy will cause the wavevectors of the two 
order parameters to be locked into equality with one another. (In this case minimization 
of 5f will also fix the relative phase (pwri — 0lti-) Since exchange anisotropy is usually not 
large, the wavevectors associated with different spin components are normally almost equal 
to one another. In that case /lock will be large enough to lock the HTI and LTI wavevectors 
to a common value. This "locked" scenario is quite common and we assume it to be the 
case here. Indeed for the systems discussed below, the data suggests that the HTI and LTI 
order parameters involve a single wavevector. 



/lock 



o\o'wvi{(l}n:if'<yir[i{—(liJi:if' + crHTi(— '?HTi)^crLTi('?LTi)^]'^< 



<?HTIi<?LTI ' 



(15) 



20 



TABLE I: General positions of the space 



CToups for NVO (top), Cmca (#64 in Ref. 



and 



4J). For Cmca the primative translation vectors are 



for TMO (bottom), Pbnm (#62 in Ref. 
ai = {a/2)i + {b/2)j; 3.2 = {a/2)i — {b/2)j; and as = ck. For Pbnm they are ai = ai; 3l2 = bj, and 
a3 = ck. 



r = {x,y,z) 


2^r = {x,y + \,z + \) 


2yT ={x,y + \,z+\) 


2^r = {x,y,z) 


It = (x, y , z) 


nixyY = {x,y+\,z + \) 


rrixzY = {x,y + ^,z+ ^) 


rriyzT = (x, y, z) 



r = 


{x,y,z) 


ruxyr = {x,y, \ - z) 


2^r = (x + i,y + i,2;) 


ma^zTC = {x + \,y + \,\ + z) 


Ir = 




2zr = {x,y, l + z) 


niyzT = - x,y + ^,z) 


'^y^ ~ (2 ~ •'^^y ~^ 2^ 2 ~ ^) 



V. MAGNETIC SYMMETRY 



A. Nontrivial Unit Cell 



There is one final refinement of our toy model which we must consider, namely the 
structure of the magnetic unit cell. In the toy model considered above, the entire spin 
structure is characterized by a wavevector and a single complex vector amplitude. However, 
more generally, the wavevector determines only how the spin wavefunction evolves from one 
unit cell to the next. Now we consider how the structure of the wavefunction within a 
unit cell is restricted by the symmetry of the crystal lattice. As a preliminary, we start by 
discussing the crystal structure of the two systems, NVO and TMO. In Table H] we list the 
general equivalent positions which define the space group operations. For NVO we choose 
the generators of the space group to be the identity, E, a two-fold rotation about the x axis, 
{x,y,z) — * {x,y,z), the x — y glide plane, {x,y,z)— > {x,y + + i), spatial inversion, 
{x,y,z) — > (x, z), and translations. For TMO the generators are taken to be E, a mirror 
z-plane, (x, y, z) — > (x, y,^ — z), the y — z ghde plane, (x, y,z) — > (| — x, y + |, z), spatial 
inversion, and translations. The magnetic sites are at the positions listed in Table |n] and 
shown in Fig. |S1 
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TABLE II: Left: Unit cell lattice positions in NVO of the Ni^+ ions carrying 5=1 (given as fractions 
of the cell dimensions a, b, and c). The "spine" sites are and the cross-tie sites are r^. Right: 
Positions of the Mn and Tb ion sites in the unit cell of TMO as fractions of the cell sides a, b, and 
c. 



rf 


(0.25, -0.13,0.25) 


r| 


(0.25, 0.13, 0.75) 


^3 


(0.75, 0.13, 0.75) 


ri 


(0.75, -0.13, 0.25) 




(0, 0, 0) 




(0.5, 0, 0.5) 





n = 1 


n = 2 


n = 3 


n = 4 


Mn 


(0,i ,0) 


(i,0,0) 




(io,i) 


Tb^ 


{x,y, I) 


(a;+i,y + i,|) 




(x + i,y + i,i) 



a) x = 0.9836 and y = 0.0810 from Ref. |45 
B. Representation Theory 



If there are n„ spins in a unit cell, then an incommensurate state will be described by 
a wavevector q and the complex- valued Fourier amplitudes iSa,r(q), where a = x,y,z and 



r 



1,2, 



, ra^ , in terms of which we write the spin wavefunctions in the form 



^,,.(R) = 5,,(q)e^^-(^+^^)+5„.(q)*e-^'i-(^+'-^) . 



(16) 



where r,- is the position of the rth site within the unit cell. For NVO r = 1, s; 2, s; 3, s; or 4, s 
for spine (s) sites and l,c, or 2,c for cross-tie (c) sites. Note that the complex amplitudes 
S{q) are defined relative to the phase, q ■ (R + r^) which would obtain if the wave were 
perfectly sinusoidal. (This convention will simplify later results.) 

We now discuss how symmetry restricts the possible values of the amplitudes Sa,T{(l) 
and how these variables are determined via diffraction experiments. The analysis of the 
symmetry of such systems in terms of their point groups is not developed. Accordingly, 
a model-independent (representation) analysis^^ is customarily invoked in such cases. // 
the magnetic ordering transition is assumed to be continuous, then the phase transition 
is signalled by an instability in the quadratic terms when the free energy is expanded in 
powers of the order parameters. In that case, the spin ordering within a unit cell will be 
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H 2 



FIG. 8: Left: The unit cell (ABCDEFGH) showing only Ni ions numbered as in Table HIl The b 
axis is perpendicular to the plane of the paper. Dashed circles represent spins in adjacent planes 
displaced from the solid symbols by ibb/2. Right: The unit cell of TbMnOs showing only the Mn 
ions (filled circles) and the Tb ions (filled triangles), numbered as in Table HIl 

determined by the critical eigenvector associated with the first eigenvalue of the quadratic 
free energy which passes through zero as the temperature is lowered. This phenomenon 
has been discussed previously in connection with wavevector selection. As we shall see, this 
analysis of the magnetic symmetry is essential to construct the allowed couplings between 
magnetic and ferroelectric order parameters. 

In order to conserve wavevector the quadratic terms in the Landau expansion associated 
with the selected wavevector q must be of the form 

f2 = J2 Ca,r;/3,r'(q)'5a,r(-q)5/3,^'(q) , (17) 

a,l3,T,T' 

where 5a,r(— q) = <Sa,T{(l)*- For /2 to be real for any choice of the complex-valued Fourier 
amplitudes, it is required that Ca,T;f3,T'{^)* = C/3^^/;a,r(q)- As in the case of phonons or 
other normal modes, the eigenvectors of this quadratic form can be labeled according to the 
irreducible representations (irreps) r„ of the paramagnetic phase which leave the wavevector 
invariant. (This group of symmetry operations is called the group of the wavevector or 
the "little group." The relevant symmetry is that of the paramagnetic phase because the 
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expansion of the free energy in powers of the order parameters is relative to this phase .^^i^^ ) 
For the orthorhombic systems NVO and TMO considered in this review, all the irreps are 
one dimensional. So essentially, the eigenvectors must be even or odd under the rotation or 
mirror (or glide) operations of the little group. 

The result of the group theoretical analysis is that one expresses the Fourier amplitudes 
<SaTi<i) in terms of symmetry adapted coordinates associated with the irrep r„. Here 

the subscript rj indicates the spin component and type of site within the unit cell. For 
instance, for NVO, (as we will see later in Table HV|) . for n = 3 or n = 4, rj ranges over five 
values, three associated with spin components on spine sites and two associated with spin 
components on cross tie sites. For n = 1 or n = 2, rj ranges over four values, three associated 
with spin components on spine sites and one associated with a cross tie spin component. If 
we write 

Sariq) = j:ui%m(;\q) , (18) 

n,r) 

then Uj^}^ is the r/th symmetry adapted basis function of irrep n in the sense that it specifies 
the a spin component for the site r in the unit cell. These basis functions are given in 
Tables HV | IV | and ED below. (Thus, for condensation of spins in NVO via irrep Table 
IV indicates that the specification of the magnetic structure requires fixing the five complex- 
valued symmetry adapted coordinates m^J^ , ) ^^s) i ^^cy ■> ^'^d rnf^.) The advantage of 
this formalism is that the quadratic free energy only couples symmetry adapted coordinates 
having the same irrep superscript n: 

h = EEFS^^^r^yi^)^ (19) 

n rirj' 

where the reality of /2 requires that Fjj^^, = F^]^. 

To deal with this free energy it is useful to introduce normal coordinates ^^"^ in terms of 
which the quadratic form is diagonal. So we set 

<H?)=E4?e^"^ (20) 

and the unitary matrix u*^"-* (which diagonalizes the quadratic free energy) is chosen so that 

/2 = T.W;^'. (21) 

n,p 
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Usually 

A(") = dkT - ~Xp{n) , (22) 

where A^"^ is a temperature-independent interaction term. As we shall see, the transforma- 
tion U to symmetry adapted coordinates is determined by the symmetries of the system, 
whereas the further transformation to normal coordinates depends on the details of the 
interactions between spins. However, the explicit form of u will not be needed for our 
analysis. 

To make this analysis more intuitive we may liken it to the problem of a particle moving in 
a spherically symmetric potential. For illustrative purposes we consider a Coulomb potential 
with a weak spherically symmetric perturbation. To solve this problem, one introduces 
symmetry adapted basis functions for the various irreps, which in this case are s functions 
ip^\ p functions ■ipn'^\ '4'n^\ ^n'^K c^c. Then the nth eigenfunction of type q {q = s, 
q = Px, Gtc.) is written in terms of the basis functions as 

^If^ = E^Vi"^^ • (23) 

m 

The basis functions ip^^ are constructed solely from symmetry arguments. The actual value 
of the transformation u to eigenf unctions depends on the details of the potential. 

We now return to the problem of incommensurate magnets. At high temperature, i. 
e. in the paramagnetic phase all the eigenvalues A|,"^ are positive and the trial free energy 
is minimized by setting ^^"^ = for all p and n, so that all the S'Q,r(R) are zero: the 
paramagnetic phase is stable against the formation of long range magnetic order. As the 
temperature is lowered, one of these eigenvalues will pass through zero [c. f. Fig. HJand Eq. 
()22|) ] and the irrep for which this happens becomes "active," so to speak. It is conceivable 
that eigenvalues A^"^ and A^? ■* of two different irreps could be degenerate. We reject the 
possibility of such an accidental degeneracy. However, if one adjusts an additional control 
parameter, such as the pressure, it is possible to reach a higher order critical point where 
two irreps simultaneously become active. A simple example of this principle arises when 
one treats a ferromagnet on a tetragonal lattice. In that case one irrep is one dimensional 
and corresponds to the ferromagnetic order lying along the four-fold crystal (c) axis and the 
other irrep is two-dimensional and corresponds to ordering in the plane perpendicular to 
the c axis. Clearly, the mean-field transition temperatures for these two distinct orderings 
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should be assumed to be different. If the anisotropy is easy- axis the ferromagnetic moment 
will lie along the c axis and if the anisotropy is easy plane the moment will be perpendicular 
to the c axis. It is possible for the anisotropy to vanish, but only by adjusting another 
thermodynamic variable, such as uniaxial stress. One therefore concludes that criticality is 
associated with a single irrep. Since the transformation to symmetry adapted coordinates 
can be determined using only symmetry considerations, the possible patterns of spin ordering 
within the unit cell are strongly restricted. 

In the usual presentation of representational analysis,— the only constraint on the sym- 
metry adapted coordinates is that they transform properly under the operations of the little 
group, Qk- For NVO the wavevector in question is along the crystal a axis, and the generators 
of Qk are 2^, a two- fold rotation about the x-axis (we often refer to a, b, and c as x, y, and z, 
respectively) and m^y, a glide operation which takes z into —z followed by a translation of 
(b + c)/2. For TMO the wavevector is along the crystal b axis and the generators of Qk are a 
two- fold screw rotation about the y-axis, 2y and the mirror plane rrixy. These operations are 
defined in Table HI It should be noted that in the case of TMO the four symmetry operations 
do not actually form a "group" because [2y]'^ picks up a phase factor and is thus not equal 
to a member of the group. This situation also occurs in connection with the application of 
group theory to the band structure of nonsymmorphic space groups . ^^i'^^ (These are space 
groups for which some pure point group operations are not space group operations.) The 
formal solution to this problem is cumbersome.— Here we may avoid these complications by 
defining the operator 2y = e^'^"'^'^2y, so that [2y]'^ translates the wave into itself. The charac- 
ter tables for these little groups are given in Table HTTl (In essence, the character table tells 
whether wavefunctions associated by r„ are even or odd under the symmetry operations 
hsted.) 

For either NVO or TMO the next step is to construct the symmetry adapted basis func- 
tions which transform according to the irreps listed in Table HTTl These allowed symmetry- 
adapted spin functions are listed in Table ITVl for NVO and in Tables IVl and IVTl for TMO. 
(In Table IVl the symmetry adapted coordinates are denoted and in Table EH they are 
denoted M^'™.) Note that each symmetry adapted coordinate appearing in these tables is 
complex-valued, each with an amplitude and (within the analysis discussed up to now) an 
independent phase. Note also that simply specifying the symmetry does not fix all the spin 
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TABLE III: Character Table for the symmetry operations O of the group^'^^ for the irreps 
r„ for incommensurate magnetic structure of (left) NVO with k = (q, 0, 0) and (right) TMO for 
k={0,q,0). 
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1 


2 


Tflxy 






C)(a) = 
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9 






ri 


1 


1 


1 


1 




Tl 


1 


1 


1 


1 


r2 


1 


1 


-1 


-1 






1 


1 


-1 


-1 


r3 


1 


-1 


-1 


1 






1 


-1 


1 


-1 


^4 


1 


-1 


1 


-1 




r4 


1 


-1 


-1 


1 



a) For an operator O we define O = e'^'^/'^O. 
b) Operators (without tildes) are defined in Tabled 

components in the unit cell. Rather it allows a choice of spin components for each set of 
symmetry-equivalent sites. For NVO one can specify the spin components of a single spine 
site and those of a single cross-tie site (unless a component is forced by symmetry to be 
zero). Having done this, the spin components of the other spine and and cross-tie sites are 
then fixed by the symmetry properties of the irrep in question. 

When the temperature is lowered further and the LTI phase is entered, then an additional 
irrep will become active via a second continuous phase transition. For NVO the new LTI 
representation (in addition to r4 already present in the HTI phase) ia^i Fi and for TMO the 
new LTI representation (in addition to F3 already present in the HTI phase) is^ F2. In an 
appendix we discuss that when two different irreps are active, their presence does not induce 
the development of a third irrep. However, had there been a further phase transition from 
the LTI phase into yet another incommensurate phase with three irreps, then the presence 
of three different irreps would induce the presence of a fourth one. 

As an illustration of how to apply the above results, we show, in Fig. typical spin 
configurations for NVO which result from the spin wavefunctions which transform according 
to irrep #4. The configuration shown in the top panel is not allowed because, as we discuss 
below, it does not respect inversion symmetry. 
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TABLE IV: Symmetry adapted basis functions UaT-r) which transform according to the irreducible 
representation r„ for the incommensurate phase associated with k = {q, 0, 0) for the Ni spine (s) 
and cross-tie (c) sites of NVO. Here Sis denotes the vector with components Sx^u, Sy^u, Sz,is, and 
similarly for the other S's. The symmetry adapted coordinates assume complex values, as 
discussed in the text. The numbering of sites is given in Table Hll The phase factors of i are chosen 
to simplify the transformation properties under spatial inversion, as is discussed below. 







U(2) 


U(3) 


u(4) 


•Sis 


/• (1) (1) • 

[imlx , mly , imlz ) 


/ (2) . (2) (2)^ 
\^sx , ITTlsy , ITlsz ) 


/• (3) (3) . (3)^ 

[imlx , mly , imlz ) 


( (4) • (4) (4)s 
V?7?.sx , ilTlsy , ?Tis2; J 


•S2s 


(■ (1) (1) ■ (l)^ 


t (2) • (2) {2)^ 

[mix ,-imly ,-mlz ) 


( ■ (3) (3) . (3)s 

[-imlx ,rnsy ,imlz ) 


^ (4) • (4) {4)x 

y-msx ^imly' ,msz ) 


•Sss 


1 ■ (1) (1) ■ 

\-imlx, mly' , -imlz ) 


1 (2) • (2) {2)^ 

[mlxl , -imsy\ mlz) 


1 ■ (3) (3) . (3)^ 
(-zmy, rrisy ,-imsz ) 


1 (4) • (4) {4)s 

( mid - imly' , rrisz ) 


•S^s 


1 ■ (1) (1) • (l)^ 
[-trrisx , -mhy , imlz ) 


1 (2) • (2) (2). 

\m\x , imly ,-mlz) 


1 ■ (3) (3) . (3)x 
[imsx ^mhy ,-imsz ) 


1 (4) • (4) (4)s 

{-mix ,-imly , mlz ) 


•Sic 


(mE\0,0) 


(mix\0,0) 


(n (3) (3)n 

{0,mcy ,mcz ) 


(n (4) (4)s 


S2c 


i-m^J, 0,0) 


(mi^\0,0) 


{0,mfy 


(0, -m^y^mcl^) 



TABLE V: Symmetry adapted basis functions Uar-r) which transform according to the irreducible 
representation r„ for the incommensurate phase associated with k = (g, 0, 0) for the Mn sites in 
TMO. The complex- valued symmetry adapted coordinates are here denoted M^. The sites are 
numbered as in Table ITTl 



Mn Site 




U(2) 


u(3) 




1 


{Mf,Mf,Mf) 


{M^,Mf,Mf) 


{M^^MfMf) 


{Mf,Mf,Mf) 


2 


{Mf,-Mf,-Mf) 


{-Mf,Mf,Mf) 


{-Mf,Mf,Mf) 


{Mf,-Mf,-Mf) 


3 


{-M^,-M^,Mf) 


{Mf,Mf,-Mf) 


{-Mf,-Mf,Mf) 


{M^,M^,-Mf) 


4 


{-Mf,Mf,-Mf) 


{-M^,Mf,-Mf) 


{M^,-Mf,Mf) 


{M^,-Mf,Mf) 



C. EfTect of Inversion Symmetry for NVO 

Up to now we only used the consequences of the symmetry of the group of transforma- 
tions which leave the wavevector invariant. However, as we have observed previously, the 
quadratic free energy of Eq. (fTTj) or of Eq. ([11^1 must be invariant under all the operations of 
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TABLE VI: Symmetry adapted basis functions UaT-r) which transform according to the irreducible 
representation r„ for the incommensurate phase associated with k = (g, 0, 0) for the Tb sites in 
TMO. The complex- valued symmetry adapted coordinates are here denoted M'^'"'. The sites are 
numbered as in Table ITTl 



Tb Site 




u(2) 


U(3) 


U(4) 


1 

2 
3 

4 


(0,0, MP) 
(0,0,-Mj'i) 

(0,0,Mj'i) 
(0,0, -MP) 


(MP,Mp,0) 
i-Mj'\Mj'\0) 
(MP,MT.i,0) 

(-mP,mP,o) 


(0,0, MP) 
(0,0, MP) 
(0,0, MP) 
(0,0, MP) 


(MP,Mp,0) 
(MP,-Mp,0) 
(MP,Mp,0) 

(mP,-mP,o) 



the paramagnetic space group. In particular, the operations for the systems we consider here 
which are not in the little group of the wavevector are those generated by spatial inversion 
I. We now discuss how inversion symmetry places restrictions on symmetry adapted coor- 
dinates. Usually when one introduces an additional symmetry, the matrix for the quadratic 
free energy becomes block diagonal. Here, the result of the additional symmetry is not to 
reduce the size of the submatrices for the quadratic free energy, but rather it places addi- 
tional constraints on the transformation matrix u of Eq. ()20|) . As we shall see, when U is 
appropriately defined, then inversion symmetry restricts all the matrix elements F^"^, to be 
real. Consequently u is real, apart from an overall phase 0„ which is associated with each 
irrep r„. 

We first analyze the situation for NVO. We need to determine the effect of X on the spin 
wavefunctions listed in Table HVl Recall that the magnetic moment is a pseudovector. That 
means that under spatial inversion the orientation of the moment is unchanged, but it is 
simply transported from its initial location at r to the transformed location, — r. Looking 
at Eq. (fT?)|l we see that spatial inversion interchanges S and S*, but since the orientation 
is unchanged, spatial inversion will not affect the spin component label a. However, spatial 
inversion does interchange Ni sublattices #1 and #3 and also #2 and #4. In other words, 
in the notation of Table llVt we have 
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FIG. 9: Two structures for spine spins in an a -c plane for NVO obtained using spin components 
from Table ITVI which transform according to irrep r4 for wavevector k = 27r/(3a), where a/2 is the 
distance between sites. The axis for the two-fold rotation 2^ is indicated. The glide plane, which 
relates spins in adjacent layers is not shown. The sites are numbered as in Table HH At the bottom 
of each panel we give the phase (in degrees) of the wave at each site. Top: for m^sx = 1, m-it^ = —i 
and the other parameters equal to zero. Bottom: for mit^ = mit"* = 1 and the other parameters 
equal to zero. As we shall see in Eq. (|32() . the order parameters mi^ and m-it^ must have the same 
phase, i. e. m^tx l^s) must be real. Only the bottom configuration satisfies this constraint. 

To see the consequences of these relations, consider the effect of the first of these two relations 
acting on the basis functions of U*^^\ for instance. This relation is 



sx 1 sy 



(25) 



which can be written as 



Arn^^x\m^i\rn<^'J) = (^i'J, miiUi'^)* • (26) 

This same analysis can be repeated for the other representations and also for the second 
relation of Eq. (|^^ . Then we see that the choices of the phase factors i in Table IIVI leads 
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to the simple result that for a = x,y, or z, and independent of representation r„ 

Jm("J(g) = m("J(g)* . (27) 

For the cross tie sites, the situation is similar except that under spatial inversion each 
sublattice is transformed into itself. Thus, we find that 

JmW(g) = mW(g)*. (28) 

So, generally for NVO we have for any representation Tn 

lm^;\q) = m^;\q)\ (29) 

where rj = s, a oi rj = c,a, where a = x,y, z. (Note that this relation does not imply inver- 
sion symmetry. If the system has inversion symmetry about the origin, then Xm^"-* (g) = 
and magnetic order can not induce ferroelectric order. Thus one can not have ferroelectric 
order if all the m's are real.) 

Now we return to the Landau expansion of Eq. to see how these relations affect the 
determination of the critical eigenvectors. First consider the situation for NVO when one 
has Eq. Recall that 

= [F!plr ■ (30) 

Now we use Eq. fl29|) to see the consequences of the invariance of /2 with respect to inversion 
symmetry, we again use the fact that /2 must be invariant under the symmetry operations 
of the paramagnetic phase,^^^*^ and spatial inversion is one such symmetry. We find that 

= EE^S4M"ng)1[xmi?)(g)] = EE^i:>?ng)^?^ 

= T.T.F:plm(;\qrmy{q) . (31) 

Thus we see that inversion invariance of f2^^ implies that F^"^^, = F^^l^. Combining this 
with Eq. ()30|) we see that these coefficients must all be real valued. This means that all 
the components of the eigenvectors of the quadratic free energy, when written in terms of 
the variables of Eq. p9|l . can be taken to be real valued. This does not mean that these 
variables must be real. Rather, since these variables are allowed to be complex, one may 
introduce an overall complex phase factor. So, the critical eigenvector, which we denote ul^^^ 
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with p = 0, has an arbitrary overall phase, in which case the amplitudes in the HTI phase 
are given as 

= auTiu^), = dHTie^'^^^'rJ'^ , (32) 

in terms of the real-valued eigenvalue components r^^^. Because we have just found that 
the matrix F*^") is real (and symmetric), the components of the eigenvector r^o are real 
valued, but, as mentioned above, since they depend on the details of the interactions, we 
do not say anything about their explicit form. Also, because we have introduced an overall 
scale factor cthti, we may require that I]^[r^'^o]^ = 1. Equation (jH^ shows that we are 
dealing with an x — y like order parameter cr = cxHTie*'^^'^^ which has an amplitude and a 
phase. (As the temperature is varied near Thti, Landau theory gives the approximate result 
c"hti ~ (^HTi — T)^/^.) In the appendix this argument (showing that the u^^j are real apart 
from an overall phase factor) is extended to include fourth order terms in the free energy. 
In analyzing experimental data, it is very helpful to realize that apart from the overall 
phase, 0HTI; all the phases of the spin amplitudes are fixed. When speaking in terms of the 
spin components, Sarig), the listing of Table HVl indicates that (for irrep ^4, for instance), 
Sx,is{(l), Sz^isiq), Sy^ic{q), and Sz,ic{(l) will all have the same phase, but (due to the factor i), 
Sy^is{q) will be out of phase with the other variables. As it happens, unless a huge number 
(several thousand) of reflections are monitored, it is impossible to use diffraction data to fix 
the relative phases with any degree of certainty. Thus, this theoretical development is useful 
to completely determine the spin structure of complicated systems such as NVO or TMO. 

We now check to see whether or not the HTI phase has a center of inversion symmetry, in 
which case, a spontaneous polarization can not be induced in this phase. We will show that 
a phase with a single representation has inversion symmetry. First of all, because we assume 
incommensurability, we can redefine the origin to be arbitrarily close to a lattice site at R, 
such that — qR is a multiple of 2n. We have already noted that Xm^ = m^*. But if is 
redefined to be zero, this implies that Xm^ = m^, which means that the spin structure has 
inversion symmetry about the redefined origin. In Fig. El we show an example of a system 
obeying Eq. p2|l which does have inversion symmetry and one having an arbitrary set of 
parameters out of Table HVl which does not satisfy Eq. (jH^ . This latter structure does not 
display inversion symmetry. Note that, as exemplified by the bottom panel of Fig. 13 it is 
possible for a structure to be noncoUinear, but to have a center of inversion symmetry. So 
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noncollinearity, in and of itself, does not guarantee having a spontaneous polarization. 

The analysis of the LTI phase is similar. Here again one can use the transformation 
properties of the order parameters under inversion to fix the phases of the spin amplitudes. 
Again, at quadratic order, one has the same result as for the HTI phase: all the LTI order 
parameters for the LTI irrep have the same phase. The analysis is extended to quartic 
order in the appendix. 



D. Effect of Inversion Symmetry for TMO 



For TMO each Mn sublattice is transformed into itself, so for the parameters of Table |^ 
we have 



(33) 



For the Tb spins, inversion transforms sublattice #1 into #3 and #2 into #4, so that for 
them one has 



(34) 



where n = 1,2 and a = x,y, or z. [M]!^'^{q) is associated with sites #1 and #4 and M^'^(g) 
is associated with sites #2 and #3.] 

The situation for TMO is slightly more complicated than it was for NVO because of the 
presence of the lower-symmetry Tb sites. In the HTI phase one can repeat the argument 
used for NVO to show that all the spin components on the Mn sites, M^, have the same 
phase. In the HTI phase the irrep for TMO was determined to be T^. Accordingly we study 
the quadratic free energy associated with a single irrep, T^. In matrix notation we have 
the quadratic free energy in terms of symmetry adapted coordinates as 



(3) 
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(35) 



where Mq, denotes the Mn amplitude M^, and Ti and T2 denotes the Tb amplitudes, Mj'^ 
and MJ'^, respectively (all for T-^). In writing this form we have used the fact that the reality 
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of /2 requires the matrix to be Hermitian. Also the matrix elements a-h are real, as can 
be shown from the arguments used previously for NVO. We now consider complex- valued 
matrix elements z„, which have no analog for NVO.— We see that the form of Eq. 
implies a contribution to /2'^'' of the form ZiM*Ti. Using inversion symmetry this term can 
also be written as Z1T2MX. Comparing this result to that of Eq. we see that Z2 = zi. 
Similarly one can show that z\ = Z3 and Zq = z^. Inversion symmetry gives no information 
on the phase of zj. Thus the matrix for f2^^ is of the form 

a b c yi yl 
h d e y2 y*2 
c e f y-i yl 
yl y*2 Vl 9 Vi 
_yi y2 ys yl h 

One can then show that any eigenvector of this matrix must be of the form 

ij = [M,, My, M„ Ti, T2] = aHTi[r., r„ r„ c, c*]e^^«" , (37) 

where is real, c can be complex, and we require the normalization J2^a + 2|cp = 1. As 
for NVO, we introduced an arbitrary overall phase 0hti- Note that M^(g) = anTi'^ae*'^"'^', 
^I'^i^l) — crHTice*'^^'^^ and M^'^(g) = ai{Tic*e'^'^^'^^ . Thus, as a result of inversion symmetry, 
the amplitudes of the two Tb sublattices are equal in magnitude, and have equal and opposite 
relative phases (from c and c*), the value of which is not fixed by symmetry. As for NVO, 
one can verify that ip is inversion invariant if (p is redefined to be zero, since then Iva = 
and Jc = (c*)* = c. 



(36) 



E. Summary 

Finally, we should emphasize that although we do not have a quantitative treatment of 
the development of magnetic long range order, we can certainly determine the magnetic 
symmetry. This information is encoded in Table IIIII For NVO, cthti is associated with 
irrep #4 and therefore is odd under a two-fold rotation about x and even with respect 
to the mirror plane taking z into —z. Likewise ctlti is associated with irrep 7^1 and is 
therefore even with respect to both these operations. For future reference we also give the 
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TABLE VII: Transformation Properties of Order Parameters for NVO (left) and TMO (right). In 
this table cr = ce*'^. and "c.c." denotes complex conjugation. Each column gives the result of 
applying the operator at the top of the column to the order parameter listed in the row. 



Order Parameter 
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TABLE VIII: Values of the symmetry adapted parameters which describe the HTI and LTI phases 
of NVO^ and TMO.'^ The uncertainty in the last digit quoted is given in parenthesis. Where 
there is no parenthesis, the entry is fixed by symmetry to be zero. For TMO the two Tb order 
parameters were assumed to have the same magnitude (as predicted by Landau theory) and the 
phase difference between the two Tb parameters in the LTI phase was found to be 1.3(3)7r. For 
NVO T = 7K is in the HTI phase and T = 5K is in the LTI phase. For TMO T = 35K is in the 
HTI phase and T = 15K is in the LTI phase. 
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b) See Ref. 



transformation properties of cthtictlti- These results are summarized in Table IVIII The 
symmetry of the LTI phase of NVO is illustrated in Fig. ^1 

For TMO the HTI order parameter cr hti is odd with respect to the mirror taking z into 
—z and is even with respect to the mirror taking x into —x. Likewise ctlti is associated with 
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irrep #2 and is even with respect to the mirror taking z into —z and is odd with respect to 
the mirror taking x into —x and these results are summarized in Table I VI II 

In Table IVlllI we give the experimentally determined values of the symmetry adapted 
parameters that describe the HTI and LTI phases of NVO and TMO. The results for NVO 



are analyzed in detail in Refs. |2^ and We will make a few brief observations here. For 
NVO the spine spins dominantly have order in the x-direction in the HTI phase from irrep 
r4 indicating that the x-axis is the easy axis. The additional order in the LTI phase due 
to irrep Fi is along the ?/-direction, as illustrated in Fig. ^| From this figure one sees that 
interactions between nearest neighboring spins in adjacent spines displaced from one another 
along either c or b are antiferro magnetic. Since the wavevectors are the same for both types 
of order, we infer that the exchange interactions are nearly isotropic. For the Mn spins in 
TMO the situation is much the same. In the HTI phase, the Mn spins dominantly have 
order in the y direction, indicating that this axis is the easy axis. In the HTI irrep (Fa) one 
sees, from Table Fv] that sites and #2 (in one basal plane) have positive y-components 
of spin and that sites #3 and #4 (in the adjacent basal plane) have negative y-components 
of spin indicating ferromagnetic in-plane interactions and antiferromagnetic out-of-plane 
interactions. In the LTI phase of TMO, the additional irrep F2 involves spins along 2;- axis 
and Table shows that for irrep #2 the components are again arranged ferromagnetically 
within basal planes but antiferromagnetically between adjacent basal planes. The fact that 
both components of spin are organized similarly suggests that the exchange interactions are 
probably nearly isotropic. 



VI. MAGNETOELECTRIC COUPLING 

A. Landau Theory with Two Order Parameters 

Now we consider the Landau expansion for the free energy, /me, of the combined magnetic 
and electric system. One might be tempted to write 

fuE = a{T - Tf)P^ + b{T - TM)\a\^ , (38) 

where a is a magnetic order parameter and, if we wish to describe a phase transition in 
which both electric and magnetic order appear simultaneously, we would set Tp = Tm- 
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FIG. 10: Schematic representation to show the symmetry of spin components for NVO. Here we 
show sections of two adjacent a-c planes. The filled circles indicate the cross-tie sites whose spin 
components are not shown, for simplicity. The x components transform according they are 
odd under the two-fold rotation 2x about the x-axis and are even under the x — y glide plane (the 
mirror plane is shown at left with a subsequent displacement A along y). (Remember that the 
magnetic moment vector involves a cross product and therefore is an axial vector: under a mirror 
operation it picks up an extra minus sign.) The y components transform according to Ti: they are 
even under both operations. Although the x and y components have different symmetry, they can 
plausibly result from nearly isotropic exchange interactions. 

There are several reasons to reject this scenario. First of all, it is never attractive to assume 
an accidental degeneracy (Tp = Tm)- This degeneracy can happen, of course, but normally 
one would have to adjust some addition control parameter (such as pressure) to reach such 
a higher order critical point. In addition, in this type of scenario magnetic and electric 
properties would not be interrelated. In NVO and TMO, in contrast, as shown in Fig, ^ 
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the electric polarization has a dramatic dependence on the applied magnetic field,— which 
such an independent scenario could not explain. 

B. Landau Theory with Two Coupled Order Parameters 

Accordingly, we turn to a formulation in which the appearance of magnetic order induces 
ferroelectric order. (The possibility that electric order induces magnetic order is not allowed 
by symmetry, by the argument in footnote 87 of Ref. 0.) So we write 

fME = axE^P^ + a{T - TM)\(y? + Vme , (39) 

where Xe^ does not approach zero and the simultaneous appearance of magnetic and electric 
order is due to the term Vme- As we have seen, the magnetic order is associated with 
a nonzero wavevector, whereas the ferroelectric order is a zero wavevector phenomenon. 
Accordingly, we are constrained to posit a magnetoelectric coupling of the form 

Vme ~ aiq)ai-q)P . (40) 

This term will do what we want: when magnetic order appears in cr{q), it will then give rise 
to a linear perturbation in P, so that P ~ XeW{(i)\'^- This argument is schematic, of course, 
and we will have to fill in the details, which must be consistent with the crystal symmetry 
of the specific systems involved. 

The minimal phenomenological model which describes the magnetic and electric behavior 
of the HTI and LTI phases is therefore written as 

/ = ^(^ ^ 7HTl)crHTl(g)crHTl(g) + ^(^ - rLTl)<TLTl(g)o-LTl(g) 

+0{\af) + ^XEP' + VME, (41) 

where 

Vme = aA,B,y(^Aiqyo-Biq)Py ■ (42) 

A,B=LTI,HTI l=x,y,z 

C. Symmetry of Magnetoelectric Coupling 

We now show that this free energy reproduces the observed phenomenology of ferroelec- 
tricity in NVO and TMO. First, of all, in the HTI phase (where (Tun = 0) Vme is of the 
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form 



Vme — X! 

7 

where bj is real. Now we use the fact that Vme has to be inversion-invariant, since it arises 
in an expansion relative to the paramagnetic phase, which is inversion-invariant We use 
2^crHTi(9) = o'hti(q')* and XP^ = —Pj to show that must vanish. Indeed, we have already 
seen, the HTI phases of NVO and TMO are inversion invariant. So for these situations 6^ 
in Eq. PHj) must be zero and no polarization can be induced in the HTI phase. 

Now we consider the situation in the LTI phase when the two order parameters chti and 
(Jlti are both nonzero. The argument which indicated that aHTi,HTi,7 = can be used to 
establish that aLTi,LTi,7 = 0. Then we write 

Vme = 5I[c7'^HTi(g)*crLTi(g) +c*cr*Ti(g)crHTi(g)]P7 . (44) 

7 

This interaction has to be inversion invariant, so we use the transformation properties of the 
order parameters under inversion to write 

Vme = IVme = - X!['^7'^HTi(g)crLTi(g)* + c*crLTi(g)crHTi(g)*]-P7 • (45) 

7 

Comparison with Eq. ()44j) indicates that must be pure imaginary: Cy = ir^, where r-y is 
real. Then 

Vme = 2^r^aHTi('?)o-LTi(g)^7sin[0HTi - 0LTi] • (46) 

7 

This result shows that to get a nonzero spontaneous polarization it is necessary that two 
order parameters be nonzero. (A similar interaction was proposed by Frohlich et alM- in their 
analysis of second harmonic generation.) Furthermore, these two order parameters must not 
have the same phase. In fact, a more detailed analysis of Landau theory shows that the 
phase difference 0hti ~ 4>lti is expected to be 7r/2. (This result comes from an analysis of 
the quartic terms. As we observed earlier, the function of the quartic terms is to enforce 
the constraint of fixed spin length. This constraint usually means that the ordering in two 
representations should be out of phase, so that when one representation gives a maximum 
of spin lengths, the other gives a minimum of spin lengths.) 

Finally, we consider how the symmetry properties constrain the spontaneous polarization. 
Look at Table IVIII There we see how the magnetic order parameters transform under the 
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at zero applied magnetic field for NVO (left) and for TMO (right). Note the different scales for P. 



various symmetry operations of the paramagnetic phase. For Vme to be an invariant, we 
see that for NVO, P must be odd under 2^ (which restricts P to be along y or z) and it 
must be even under rrixy (which restricts P to be along x or y). Thus, symmetry restricts 
P to be only along y. This is exactly what experiment shows. For TMO, P must be odd 
under both 2y and rrixy Thus, symmetry restricts P to lie along z ai H = 0, as is observed 
in experiment. (At higher magnetic fields the magnetic symmetry must change to explain 
why the polarization switches from the z-axis to the x-axis.) Furthermore, the temperature 
variation of P, shown in Fig. ^2 looks very much like that for an order parameter. But that 
is to be expected because if we minimize the total free energy with respect to P, using Eq. 
(I46p . we see that the spontaneous polarization is given as 



-P7 ~ XeO'htiO'lti 



(47) 



When the LTI phase is entered, chti is already well developed and is therefore essentially 
independent of temperature. Thus we expect that crudely ~ ctlti- Indeed, although we 
have not undertaken a quantitative analysis, the experimental curves of P versus T look 
quite similar to those for an order parameter. 



Finally, for TMO for a large magnetic field along a (see Fig. E)) or along b (see Ref. 



4lh . there is a change of orientation of the spontaneous polarization P to lie along a. Since 
there seems to be no analogous phase transition within the HTI phase, we attribute this 
reorientation to a change in the LTI spin state. Instead of the additional irrep of the LTI 
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phase being r2 (as it is at low field), we infer that the new LTI irrep is Fi, since this 
combination of irreps is consistent with having P along a. Furthermore, if we assume 
that the exchange coupling is isotropic, then we would expect that Fi ordering would be 
ferromagnetic within basal planes and antiferro magnetic between planes. From Table IVI this 
constraint can only be satisfied if the ordering involves the x-component of spin. So, from 
the polarization data we speculate that the Mn spin structure (which at low field is in the 
b-c plane) is rotated, at high field, into the a-b plane. 

D. Broken Symmetry 

We should also mention some considerations concerning broken symmetry for NVO. 
(Clearly a similar discussion applies to other similar systems.) Since both transitions in- 
volving the HTI phase involve broken symmetry we assert the following. At the level of the 
present analysis when the temperature is reduced to enter the HTI phase, the modulated 
order appears with an arbitrary phase 0hti- Of course, if this state is truly incommensurate, 
then the phase will remain arbitrary. Normally, however, we would expect some perturba- 
tion to break this symmetry and this continuous symmetry should be removed. However, we 
do expect a degeneracy with respect to the time-reversed version of the ordered HTI phase. 
In that case upon performing many runs of the same experiment, both versions of the HTI 
phase should occur with equal probability. 

One can make much the same observation about the HTI-^LTI phase transition. Here 
one has the additional broken symmetry associated with the irrep Fi. When the temper- 
ature is reduced to enter the LTI phase, the system will have two symmetry-equivalent 
states into which it can condense. As with the usual magnetic phase transitions, one can 
(in principle) select between these two phases by applying a suitably spatially modulated 
magnetic field. Such an experiment does not seem currently feasible (because modulation of 
an applied field on an atomic scale is difficult to produce). However, because the magnetic 
order parameters are coupled to the ferroelectric moment, one can select between the two 
symmetry equivalent possibilities for the LTI order parameter by applying a small electric 
field. A interesting experiment suggests itself: compare the magnetic state as determined 
by, say, neutron diffraction for the two cases of a small applied electric field in the positive 
and negative b directions. According to the magnetoelectric trilinear coupling, application 



41 



of such an electric field should select the sign of the product cthtiClti- In this context we 
remark that measurement of the spontaneous polarization P (as in Fig. is made by 
preparing the sample in a small symmetry-breaking electric field Eg, which is removed once 
P becomes nonzero. The ferroelectric order is confirmed by verifying that P changes sign 
when the sign of Eq is changed. 

VII. MICROSCOPICS 

Since the spontaneous polarization P must result from a spontaneous condensation of an 
optical phonon having a dipole moment, we are led to study the symmetry of the phonon 
excitations at zero wavevector. Neglecting nonzero wavevectors, we write the a-component 
of the displacement of the rth ion in the unit cell at R as 



where (r) is the normalized form factor of the zth generalized displacement whose ampli- 
tude is Qi. A comprehensive analysis is given elsewhere,"^ but here we confine our attention 
to generalized displacements in NVO which transform appropriately (like a displacement 
along b) to explain experiments. Such a generalized displacement Qi must be invariant 
under the operations (see Table II) E, cr^, (Tz and 2y and change sign under ay, 2^, 2^, and 
I. There are 12 such generalized displacements of the 13 ions in the primative unit cell. Six 
of these are the uniform displacements along b of all crystallographically equivalent sites 
of a given type, viz. Ni spine sites, Ni cross-tie sites, V sites, Oi, O2, and O3 sites, and 
these uniform displacements, denoted Qi, Q2 ... Qe, give rise to a dipole moment along the 
b axis. Other generalized displacements involve, perhaps surprisingly, oppositely oriented 
displacements along the a or c axis within a group of crystallographically equivalent sites. 
We illustrate one of these (Q7 involving Ni cross-tie sites) along with Q2 in Fig. ^1 Since 
Q7 has the same symmetry as Qi ... Qq, it must couple to these modes. One can easily 
visualize this by imagining the ions to act like hard spheres. In that case, as the cross-tie 
ions approach spine sites #1 and #4, they cause these site (which initially were at negative 
y) to move to more negative y. Similarly, as the cross-tie sites move away from sites #2 
and #3, these ions have more room to move closer to y = 0. In other words, the opposing 
motion of the cross-tie sites in mode Qj along the c axis interacts with the uniform motion 




E 




(48) 
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y=-0.13b 



FIG. 12: Pattern of two generalized displacements, Q2 (left) and Qj (right), which transform under 
the symmetry operations of the crystal like a displacement along b. In the left panel "+" indicates 
motion along the positive b axis. In Q2 all the cross-tie sites move in parallel along the b axis and 
therefore this motion induces a dipole moment. As discussed in the text, the nonuniform motion 
of the cross-tie sites in the generalized displacement Q7 induces uniform motion of the spine site 
in the b direction which in turn produces a dipole moment. 



in mode Qi of the spine sites along b. In summary, the elastic free energy as a function of 
displacements can be written as 



fm}) = \f: V,,Q,Q, . (49) 

At the time of this writing no calculation or neutron experiments to determine Vij have 
appeared. Instead we have recourse to a very crude toy model, obtained by setting 

fm}) = e^^^^dQ?, (50) 

where Mj is the mass of ions in mode Qi and a;D is the Debye frequency, characteristic of 
phonons. 

We now consider the effect of a generalized displacement Qi on the exchange interaction 
between nearest neighbors in the same spine. Then for spins numbered 1 and 4 in a unit 
cell we have the exchange interaction as a function of displacement as 

Hi4(g.) = 7^14(0) -fQ,j:^-(i)^^%r^^M4) • (51) 
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The existence of a mirror plane perpendicularly bisecting the 1-4 bond (mbc) implies that 



(52) 



which is 
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(53) 



where we used mi,cQi = Qi, "rribcSix = Six, m^cSiy = -S^y, and mbc^i^ = -S^z (the spin 
is a pseudovector) . Thus the derivative must have the form 



dQi 
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(54) 



where J^^ is the symmetric exchange tensor and D is the Dzialoshinskii-Moriya vector, 
which specifies the antisymmetric component of the exchange tensor. 

We determine the other similar interactions in the unit cell using the appropriate sym- 
metry operations. If 2y is a rotation about an axis parallel to b and which passes through 
site #4, then 



H(4,l';y) = 2,7^(1, 4; y) 
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where we used 2yQi = Qi, and site 7^4' is one unit cell to the right of site #4 in Fig. [T21 
Also, if 2x is a rotation about the a axis, then 



n{2,3-y) = 2,n{iA\y) 
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where we used 2xQi = —Qi and 

n{3,2';y) = 2,H(2,3;y) 
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dQi 
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(57) 
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where site #2' is one unit cell to the right of site #3 in Fig. IT^ 

When we consider Eq. ISUj) and neglect fluctuations, the spin-phonon interactions lead to 
the result 



(58) 

/II r ■ 

where ( ) indicates a thermal average and [n, m) are summed over the 4 nearest-neighbor 
spine-spine interactions in a unit cell. Assuming the spins are characterized by spine spin 
components scaled by a for irrep 4 and by b for irrep 7^1, we write the spin components as 
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Sz{xi) 

Sz{x2) 
Sz{x3) 

SJxA 



(a, + zfe.le'""! + (a! - ib^e 



iqxi 



{a, - i6,)e^''^i + {al + ibl)e-"'''^ 
(a, - «6,)e^^"^ + {al + z6:)e-'«"^ 
(a, + z6,)e'«^'^ + {al - ibDe-"^'^' . 

Using these evaluations one can carry out the sum over (ra, m) in Eq. ()58|1 to get 

(g,) = 16(M,a;^)-i[F/^)sin(ga/2) +Ff)cos(ga/2)] , 

where 

f/""-* = 5>[a*6^ + albx]dDy/dQi + E TTa^iaabl.jdJaa/dQi (63) 

a 
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and 

F^^'^ = Q[a,b; + bXyWyJdQi + Q[a,b; + b,a;]dDjdQi , (64) 

where —ir^ = TCy = tc^ = ^- Note that these terms require the presence of both order 
parameters a and b and hence they can only be nonzero in the LTI phase. Also these 
terms are only nonzero if a and b have different phases. For displacements which could 
give rise to a spontaneous polarization along the a or c axes, the sum over (n, m) in Eq. 
fl58|) gives zero.— These conclusions agree with the result found using Landau theory. This 
magnon-phonon coupling also contributes to the temperature dependence of the wavevector 

^ 22.23 

To get an order-of-magnitude estimate of the various quantities, we consider the effect of 
the motion of the oxygen ions, which are the lightest atoms and therefore have the largest 
displacements. Crudely speaking, the dipole moment, Pq of the generalized displacement Q 
is given by Pq = qg J2ieQ Ui, where gg is the charge of the ions of Q and Ui is the displacement 
of ion i in Q. More accurately, qq should be replaced by an effective charge qg which takes 
account of the electrical relaxation which occurs as the ions move. (This is analogous to 
the discussion given at the end of the preceding paragraph.) Thus even Qj will develop a 
(probably small) dipole moment in the b direction. However, for simplicity we set 

(P) = qi{Q)/vnc = 2e(Q)Kc , (65) 

where ~ 275 x 10~'^°m'^ is the volume of the unit cell. More accurately, (Q) should 
be replaced by {Q)\/n, where n is the number of ions involved in the mode generalized 
displacement Q. (So n = 4 orn = 6.) So, in meters, (Q) ^ (275 x 10"^°) (5 x 10"^)/(3.2 x 
lO^^^y^), where we took P = 5 x 10~'^C/m^ as a typical value. Thus we estimate the ionic 
displacement to be of order (Q) ~ O.OOlA. (Actually, neutron diffraction indicates that the 
ionic displacement ought to be at most O.OOlA.—) If dJ/dQ represents a typical value for 
dM^p/dQ, then 

{hc)"^ dJ , . 

Working in A and eV and taking (Qi) = O.OOlA, hun ~ 0.05eV, M^c^ ^ 1.6 x 10^°eV, and 
he ~ 2000eVA, we find that this mechanism requires that 

dJ 

— ~ O.OleV/A • (67) 
oQ 
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This seems to be a plausible value. Obviously a first principles calculation of dMap/dQi 
would be of interest to make this analysis more concrete. 

VIII. SUMMARY AND OUTLOOK 

The development of multiferroic materials having very large magnetoelectric couplings 
offers the possibility of designing new types of devices which exploit the coupling between 
magnetic and ferroelectric order. Furthermore, investigating the nature of the coupling 
between magnetic and ferroelectric order parameters in these compounds may be important 
in understanding other systems displaying significant interactions between different types of 
long range order. We will briefly summarize the main results of the model we have presented 
coupling ferroelectricity with incommensurate magnetic order, and then discuss what this 
implies for future research on magnetoelectric multiferroics. 

A. Summary of this Review 

As many of the recently identified materials exhibiting simultaneous magnetic and ferro- 
electric order are incommensurate magnets, we have focused on these systems. We discussed 
a toy model for incommensurate magnetism. In this model, we saw that, under the assump- 
tion that the magnetic anisotropy is not too large, the magnetic system will undergo a para- 
magnetic to a longitudinally ordered incommensurate phase we refer to as the HTI phase. 
On further lowering the temperature, the is another transition to a distinct incommensurate 
phase with additional transverse spin ordering, which we call the LTI phase. 

Because knowing the detailed symmetry of the incommensurate magnetic structure is 
crucial for determining whether or not ferroelectric order is allowed, we considered the 
extension of this toy model to systems with non-trivial unit cells. We addressed this problem 
by expressing the spin order parameters in terms of irreducible representations consistent 
with the symmetry restrictions of the unit cell. The central observation for understanding the 
magnetoelectric coupling is that the free energy must be invariant under all symmetries of the 
paramagnetic phase, and in particular, if the paramagnetic crystal has inversion symmetry, it 
must be invariant under spatial inversion. This requirement was used to determine whether 
or not a particular incommensurate magnetic structure allowed the possibility of ferroelectric 
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order. Using this approach, we are able to quahtatively explain the multiferroic behavior 
of both NVO and TMO, including the absence of ferroelectric order in the HTI phase, 
the development of ferroelectricity in the LTI phase, and the qualitative features of the 
spontaneous polarization (direction and temperature dependence). 

It is worth noting that the magnetoelectric coupling we have described here does not 
reduce to the analogous coupling which can occur in a ferromagnet or in an antiferromagnet. 
As remarked in the review of Smolenskii and Chupis,— such a trilinear coupling can not exist 
in structures which (like NVO or TMO) have inversion symmetry in the paramagnetic phase. 
Indeed, the mechanism we invoke requires that g 7^ 0, as one can see from Eq. ()4fi|l . (If q = 0, 
then the order parameters are real, 0hti = 0lti = and Vme = 0.) In that review they also 
mention a coupling which involves gradients of the magnetic order parameter. That type of 
coupling may be related to that used here, although in our case the symmetry properties of 
the unit cell play a crucial role which can not be replaced by a continuum vector field. 

We also showed that the microscopic symmetry of the derivative of the exchange tensor 
with respect to ionic displacement leads to results in complete agreement with the symmetry 
arguments based on the Landau expansion. This symmetry will have to be respected by any 
truly microscopic theory of magnetoferroelectrics. 

B. Outlook for Device Applications 

The success of the theory described in this review suggests that it may be valuable both in 
understanding the origins of multiferroic behavior in presently identified systems, and well as 
in guiding the search for new multiferroic compounds having desirable materials properties. 
We briefly discussed the technological drivers motivating the search for magnetoelectric 
materials, by illustrating the types of devices that might be possible using multiferroics. 
However, one of the difficulties that must be resolved before these systems could possibly be 
incorporated into fabricating next generation magnetoelectric devices is that the very low 
transition temperatures into the ferroelectric LTI phase (6.4 K for NVO and ~ 27 K for 
TMO) make these materials unsuitable for many applications. We discuss in the following 
some of the general ideas extracted from our model which may help guide the search for new 
multiferroics having higher transition temperatures and larger magnetoelectric couplings. 

Extending the search for multiferroics from simple ferromagnets to systems with incom- 
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mensurate magnetic order is an important first step in finding materials which have a room 
temperature transition into a phase with simultaneously appearing magnetic and ferroelec- 
tric order. Insulating ferromagnets tend to have very low transition temperatures, but 
many incommensurate magnets have ordering temperatures well above room temperature.— 
In fact, the incommensurate magnetic structure associated with ferroelectric order in one 
recently identified multiferroio^ persists up to T = 320 K. Our results suggest that insu- 
lating incommensurate magnets with high magnetic ordering temperatures may be prime 
candidates in the search for strongly coupled magnetoelectric multiferroics at room temper- 
ature. 

In the next subsection we propose experiments which might show that in systems such as 
those studied here magnetic phase boundaries can be sensitive to the applied electric field. 

C. Experimental Outlook 

It is suggestive that for the two systems we have considered in detail, the development of 
a single order parameter can not induce a spontaneous polarization. A requirement for such 
a result to be general is that in the paramagnetic phase the system should have a center of 
inversion symmetry. In that case, it is assured that in the paramagnetic phase there is zero 
spontaneous polarization. We speculate that for such a system, the development of incom- 
mensurate magnetic order via a continuous phase transition does not allow a spontaneous 
polarization to be induced by the magnetic ordering. That is, the development of a single 
magnetic order parameter transforming under one specific irrep cannot break the spatial 
inversion symmetry and lead to ferroelectric order. In order for incommensurate magnetic 
order to break the spatial inversion symmetry (and potentially induce the development of 
ferroelectricity) we postulate that it is necessary to have two distinct magnetic order param- 
eters. (A similar proposal in connection with second harmonic generation has been made 
by Frohlich et alM-) 

This speculation raises the more general question concerning a system which in the dis- 
ordered phase has neither magnetic nor ferroelectric order. In particular, consider systems 
which lack inversion symmetry, but whose rotational symmetry elements preclude a nonzero 
vector order parameter. We give two families of such crystal structures. The first is that 
of the point group D2 (orthorhombic space groups #16-^24 in Ref. |44 ) and the second is 
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that of point group T (cubic space groups #195 - #199 in Ref. 144 ) . In the para phase, 
these systems have no magnetic long-range order and, because these crystal structures do 
not allow vector ordering, they do not display ferroelectric order. When such a system de- 
velops long-range incommensurate order with a wavevector along one of the crystallographic 
directions, then only rotations about this direction remain symmetries and a spontaneous 
polarization along the direction of the wavevector is permitted, at least in principle. So in 
this case, it would seem that two order parameters would not be necessary for magnetic long 
range order to induce ferroelectric order. Accordingly we suggest that it would be interesting 
to find ferroelectric incommensurate magnets having one of these crystal structures. 

Another experimental program which this study suggests concerns the phase diagram of 
these systems in the T-E plane, where E is the uniform applied field. Since only the LTI 
phase of NVO or TMO has a spontaneous polarization, this phase is favored (relative to 
the HTI or AF phases) in the presence of an electric field. So we propose the schematic 
phase diagram shown in Fig. Uni We have not specified the scale of the horizontal axis 
in this schematic figure, but at least for the LTI-AF transition in NVO we can estimate 
how thin a film would have to be to produce a 5% shift in the transition temperature for 
an applied voltage of 5 V. The analog of the Clausius-Clapeyron equation for the LTI-AF 
phase boundary in the T-E plane is 

dT/dE = -ViPi^r, - Paf)/(5lti - Saf) • (68) 

Now take V to be the volume per Ni ion. The volume of the conventional unit cell is 
V = abc, so V = abc/12, because there are 12 Ni's per conventional unit cell, a = 5.9, 
b = 11.4, and c = 8.2, all in lO^^^m. For one Ni ion the total entropy change from zero to 
infinite temperature is fcln3. Guided by specific heat measurements^! we set 

(5+ - 5_) = 0.01A;ln3 (69) 

and take P = 5 x lO^^C/m^ as a typical value. Then we find 

dT/dE = -[(540 x 10-^°mVl2])(5 x 10-^C/m^)/[1.4 x lO'^^J/K] 

= 1.5 X 10"^Km/V . (70) 

To see roughly what this means, set dT = 0.05Tc = 0.2K and dE = 5V/t, where t is the 
thickness of the sample. This gives 0.2t x 10~^, or t ^ l/x. 
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FIG. 13: Schematic phase diagram for A system hke NVO in the T-E plane. (We do not give the 
scale of the £^-axis.) 

Finally, we emphasize that it would be desirable to determine the magnetic structure 
of TMO for high magnetic fields along the a or b direction to test whether the magnetic 
structure proposed below Eq. (jT7|) is realized. 



D. Theoretical Outlook 

It is clear that the next step for theorists is to construct a fully microscopic theory to 
explain the phenomenological trilinear interaction highlighted in this review. Here we in- 
dicated how the dependence of the exchange interaction on ionic displacements gives rise 
to the symmetries expected from Landau theory. What is clearly missing is a microscopic 
calculation of the exchange constants. This sort of calculation as a function of bond angles 
has been pursued for Cu-O-Cu bonds.— However, what is needed here is the more compli- 
cated calculation for Ni-O-Ni bonds and furthermore, it would seem that this is going to 
require some sort of calculation based on the local density approximation to determine the 
dependence of the exchange tensors on ionic displacements. Calculations of this type are 
being carried out. 
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E. Closure 

Magnetoelectric materials have been investigated in depth for the last forty years. How- 
ever, we are still identifying completely new classes of ferroelectric materials, showing new 
and different types of magnetoelectric couplings. With the increasing interest in materi- 
als which couple electric and magnetic properties, we expect that these magnetoelectric 
multiferroics will continue to be an active area of research. 
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APPENDIX A: QUARTIC TERMS 
1. HTI Phase of NVO 

Since irrep #4 can not induce any other irrep, the free energy only involves order param- 
eters of that irrep. (So we will omit the superscript n = A which labels that irrep.) Then, 
correct to quartic order we write the free energy associated with irrep #4 as 

n = n^^^ + H^"^ + . . . , (Ai) 

where 

tt' 

and 

7^(4) = (1/4) {qymr^{qymr.,{q)mr^{q) , (A3) 

where m^=i(g) = m^s{q), mr=2{q) = mys{q), m^=3(g) = m^s{q), rnr=i{q) = my^q), and 
m^=5(g) = rrizc, and the spin components at all the sites in the unit cell are given in Table HVI 
in terms of these variables. Hermiticity implies that Vr,T' = v*,^^ and Wri,T2,r3,r4 = '^Ta,T4,Ti,T2 
and w can be taken to be symmetric under interchange of ri and T2 and of T3 and r4. As 
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discussed in the text, using inversion symmetry one can show that all the matrix elements 
of V and w are real. 

Now we transform to normal modes: 



=T.^rAp^''"' , (A4) 



^p' = E rr'p'mAq)e-'^^' ^ ip'e~"^^' , (A5) 

t' 

where p = 0, 1,2,3,4 labels the normal mode, the r's are real, and the critical mode (po) 
has an amplitude which heretofore we called cthti- The quartic Hamiltonian is 

'^4 — T '^ri,T2,T3,T4l^Ti,pi^T2,P2 
n,T2,T3,T4 pi,p2,P3,P4 

xr r f f f f g»(</'P3+9^P4-<?^pi-'^P2) (AQ) 

This quartic term will involve contributions proportional to C,o = c^htI' where p ranges from 
zero to four. If we were to omit the quartic terms with p = 3, then the minimum of the trial 
free energy would be realized for (Thti 7^ 0, but with the other ^p's being zero. Therefore, 
the most important term to consider is the one cubic in cthti, which is 
1 ^ 

S'Hi = — E E ^'^l>f"2.T"3,T4^ri,0^T2,0^T3,0'"r4,pO^HTI^ + C.C. 

T1,T2,T3,T4 p=l 

3 

= J2 ^pO^HTI COs[0p - 0HTl] , (A7) 
p=l 

where Ap is real. The quadratic terms for the noncritical variables can be written as 

^ p=l ^ p=l 

where Xp is the susceptibility of the pth mode. Then, after minimization with respect to the 
noncritical variables we see that cos[0p — 0hti] = ±1 (so that (Tlti and ctlti are in phase^S) 
and 

^p = XpApaliTi , P > . (A9) 

Thus the effect of the quartic terms is to induce nonzero values for the noncritical normal 
modes and thereby slightly change the components of the critical eigenvector, but the quartic 
terms do not change the fact that all the order parameters belong to irrep 7^4 and that 
they all have the same relative phase. 
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2. LTI Phase of NVO 



Now we consider the LTI phase, where we have two irrep simultaneously present. There 
are various types of quartic terms. First, consider those quartic terms which only involve a 
single irrep. We can apply the analysis of the HTI phase, to state that such terms do not 
modify the conclusion that all the symmetry adapted coordinates of irrep 7^4 have the same 
phase, 04, and all the symmetry adapted coordinates of irrep #1 have the same phase 0i. 

Next consider the more general quartic terms which involve both irreps. Terms of the 
type [m^'^^]*[m^^^]*m^^^m^^^ are independent of the phases and therefore after minimization 
of the trial free energy these terms do not modify the phases. There are no terms which 
involve three order parameters of one irrep and one order parameter of the other irrep. So 
the only terms which might affect the phases are terms of the form 
and its complex conjugate. So we consider quartic terms of the form 

Hermiticity requires that w[f!}^^,j^J^^^ = [w^l^^^^!^^^^^* . Then inversion symmetry indicates 
that the w coefficients are real. Thus these quartic terms give 

= Acos[2(04-0i)] E E -SL-S.-S3-2UJ?eJ?ee ' (All) 

where the r's are the real- valued transformation coefficients determined in quadratic order. 
All the quantities in F4 are real. So F4 is minimized by either setting cos[2(04— = ±1. An 
explicit calculation for the actual experimentally determined values of the order parameters 
indicated that the correct choice of sign is the negative sign, and therefore that the two 
irreps are out of phase with one another. This conclusion agrees with the intuitive argument 
based on the idea that quartic terms tend to enforce the fixed spin length constraint. When 
the coordinates of one irrep are maximal, then those of the other irrep should be minimal. 
Thus we conclude that | sin(0HTi — 0lti)| = 1 in Eq. 

APPENDIX B: DO TWO IRREPS INDUCE A THIRD ONE? 

When two irreps, r^(g) and Ty{q) are simultaneously present (as happens in the LTI 
phase), one might ask whether their combination could then induce a third representation. 
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Taiq), all of which are assumed to be associated with the selected wavevector q. Since 
ry{q)ry(—q) is unity, it is equivalent to ask whether for some k, products hke T^iq)''''^ <8) 
ry{—q)''^Ta{q) or Ty{q)''~^<SiTx{—q)''(^Ta{q) transform like unity. (The form of this product 
is dictated by wavevector conservation. In this connection we neglect the possible effects of 
Umklapp terms.) If one of these products satisfies this condition, then the existing order 
parameters can give rise to a linear field acting on Taiq), thereby inducing a nonzero value 
for this representation. By explicit enumeration of the various cases one can verify that the 
condition to induce a third irrep can not be satisfied. If, hypothetically, there existed a third 
phase transition in which a third irrep condensed, then the presence of these three irreps 
would induce the fourth irrep. 
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